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Abstract 

We have constructed a spherically symmetric structure model within an otherwise FRW uni- 
verse filled with perfect fluid and studied its quasi-local characteristics. This is done by using 
the Lemaitre solution of the Einstein equations and by suggesting an algorithm to integrate it 
numerically. The result shows how a dynamical black hole evolves within the FRW universe, 
its decoupling from the expanding parts of the model, the structure of its space-like apparent 
horizon, the limiting case of the dynamical horizon tending to a slowly evolving horizon, and 



X 

, the decreasing mass in-fall to the black hole. The quasi-local features of this cosmological black 

hole may not be inferred from the Newtonian approximation although the gravity outside the 
structure is very weak. 
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I. INTRODUCTION 



The term cosmological black hole is used to describe a collapsing structure within an 
otherwise expanding universe. Since the early beginning of the discovery of the expansion 



of the universe people have 
a cosmological background ( 



Deen looking 



l|, see also {2]). Such an overdense region, whcih may be 



br models describing an overdense region in 



expanding at the beginning, will finally decouple from the expanding regions of the 
universe and collapses to a dynamical black hole different from the familiar Schwarzschild 
one despite the very weak gravitational field outside the structure. Therefore, such 
cosmological structures, if based on exact solutions of general relativity and not produced 
by a cut-a„d-paste technology, are very interesting laboratories to stndy not only^eneral 
relativistic structures, their quasi- local features such as mass and horizons [3], but 
also the validity of the Newtonian approximation in the presence of the very weak 
gravitational fields related to non-local concepts within general relativity. After all, the 
universe is evolving and one needs to have a dynamical model for a black hole to be 
compared with the familiar results in the literature on black holes within a stationary 
and asymptotically flat space-time where global concepts such as event horizon are 
not defined. The need for a local definition of black holes and their horizons has led 



horizon 



a 



us to concepts such as Hayward's trapping horizon 5|], isolated horizon joj, Ashtekar 
and Krishnan's dynamical horizon (DH) [7j, and Booth and Fairhurst's slowly evolving 



Now, a widely used metric to describe the gravitational collapse of a spherically sym- 



metric dust cloud is the so-called Tolman-Bondi-Lemaitre(LTB) metric {q]. It was pointed 

ack holes. These models may be ex- 



Our interest is now using these 



out in j^, Q that the model admits cosmological b 
tended to a perfect fluid with a non zero pressure 
cosmological solutions to construct dynamical black holes within a FRW expanding uni- 
verse, and study their characteristics. To do this, we have to avoid any cut-and-paste 
method of finding the solution. Section II is an introduction to these inhomogeneous pre- 
fect fluid cosmological models and how numerically integrate the field equations leading 
to a dynamical structure within an otherwise FRW universe. In section III the result of 
numerical integration is reported expressing the main characteristics of our models assum- 
ing different pressure profiles. We will then discuss the result in section IV. Throughout 
the paper we assume SnG = c = 1. 



II. GENERAL SPHERICALLY SYMMETRIC SOLUTION 



Consider a general inhomogeneous spherically symmetric spacetime filled with a perfect 
fiuid and a metric expressed in the comoving coordinates, = {t, r, 6, </>): 

ds^ = -e^'' dt^ + e"" dr^ + RUn\ (1) 

where a = a(t,r), X = X(t,r) are functions to be determined, R = R{t,r) is the phys- 
ical radius, and dfl"^ = dO"^ + sin^ dcj)"^ is the metric of the unit 2-sphere. The energy 
momentum tensor of the perfect fiuid is given by 

T^- = (p + p)M^M- + ^A«-p, (2) 

where p = p(t, r) is the mass-energy density, p = p{t, r) is the pressure, and = 
{e~" , 0, 0, 0) is the perfect fiuid four-velocity. 



A. Field Equations 

In addition to the Einstein field equations, G^^ = nT^^ — g^'^A, we will use the 
conservation equations in the form 



VJT^^ = A + - — '-- + — = 



{p + p) 



(p + p) R 



(p + p) 



V„T''^ = a' 



p 



p + p 



0, 



(3) 



(4) 



where the dot means the derivative with respect to t, and the prime means the derivative 
with respect to r. The Einstein equations lead finally to the two following equations 



d_ 
dr 



and 



d_ 
di 



R + RR^e-^" - RR'^e-^ - -AR^ 



1 



KpR^R', 



-KpR^R. 



R + RR^e 

The term in the brackets is related to the Misner-Sharp mass, M, defined by 

R 3 



Eqs ([5]) and (|6]) may now be written as 



np 



2M' 

Wr 
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2M 



(5) 



(6) 



(7) 



We may write Eq ([7]) in the form of an evolution equation of the modeh 



where 

f{t,r) = R"e-'-l (10) 

is the curvature term, or twice the total energy of test particles at r (analogous to /(r) 
in the LTB model). Note that R{t,r) can not be directly obtained from this equation 
because of the unknown functions A, a, and M. 

The metric functions gu and Qrr may be obtained by integrating ([3]) and (j4]): 



p' dr 



a - c{t) - ^|.^co„. - ao - J^^ j^f^^ dp (11) 



A = Ao(r) - 2 / - 41n f ^ ) \t=consu (12) 



and 

\ = \^(r'] - 9 1 , 

(p + pip)) \Ro 

where {t) and Xo{r) are arbitrary functions of integration. In the case of c{t) it is easily 
seen that requiring our coordinates to lead to the LTB synchronous ones for p = leads 
to c(t) = 0. We notice also that according to (ITOl) and the LTB coordinate conditions, the 
choice of Ao(r) is equivalent to the choice of f{tQ,r) = fo{r). One may prefer to choose 
/o(r) and then calculate Ao(r) from e^" = -Rq /(I + fo). 

We have therefore 5 unknowns p, p, a, A, and R, four dynamical equations p, M, A, and R, 
in addition to an equation of state p = p{p), and the definition of the mass M ([7]). This 
defines a numerical algorithm to find solutions for the dynamics of our spherical structure 
after assuming some initial conditions. 



B. Construction of the Lemaitre Model 



To generate a Lemaitre model in the general case we need a numerical procedure. 
We first specify the arbitrary initial functions, -Ro('^) = R(to,r), Po(^) = p{to,r), Ao(r) = 
X(to, r) , cro(r) = cr(to, f), and the equation of state p{p). We then integrate the dynamical 
equations for constant t or r in the following order: 

1. Choose an initial time to, and specify R{to,r) = Ro{r); Rq is then also known form 
the derivative of (fT5]l with respect to r at t = to! 

2. Specify M{to,r) = Mo(r) at t = to; 



3. Once Mo(r) and Ro{r) are specified, po{r) may be determined from (jS]); 

4. Select an equation of state, p = p{p); 

5. Choose A(to,''"), or choose first f{tQ,r) = /o(r) and then calculate Ao(r) from e^'^ = 

^[,7(1 + /o); 

6. (T(to,r) can then be obtained by integrating (fTTj) along t = to- 

We have now specified how to determine all the needed initial functions at the time to- 
Therefore, their r— derivatives are also known. The time evolution of the metric functions 
along the worldlines of constant r may then be calculated. This is done in the following 
way: 

1. From equation f|T3|) we obtain R: 



2. Eq (JED then give us M: 



i? = ±e^^/ — + / + ^ ; (13) 



M = ; (14) 



3. Combining Eq ([3]) and Gqi allow us to eliminate a'. Then substituting A from 
we arrives at p: 

(15) 



R 

p = -p' ^~{p + p) 



R' 2R 



4. From the equation of state p = p{p) we obtain p: 



p=—p; (16) 
dp 

5. Eq (jlj) may be combined with (fT5|) to give 

6. Using the initial values for t = t^, we then are in a position to solve the above 
5 differential equations numerically to obtain R(t,r), M{t,r), p(t,r), p(t,r) and 
A(t, r) for every t and r. Note that in each step the spatial derivatives R', p' and 
M' needs to be determined; 

7. Finally, oit^r^ is obtained from Eq ( ITTi) by integrating along constant t. 

Notice that we have chosen 4 initial functions, -Ro('^), Po('^), '^o(^) and o"o(t), as well as 
the equation of state ]?(p). 



III. INITIAL PROFILES AND THE RESULTS 



We are now ready to specify the equation of state and integrate the model to see its 
characteristics. We consider two types of fluid which helps us to have a comparative 
discussion of the result: a perfect fluid with a constant of state function, p — wp, and 
the more general case with the equation of state p — ws{r)p matching our needs of a 
structure with pressure inside and a pressure-less matter dominated universe far from 
the structure. We may then choose the function s{r) in the way that the pressure 
become zero at infinity, r >> 1. A suitable choice is s(r) = e ""o . This is a more 
realistic model to describe a black hole collapse in the FRW universe and sec the ef- 
fect of the inside pressure while the universe outside is matter dominated with no pressure. 

The model we envisage starts from a small inhomogeneity within a FRW universe. 
The density profile should be such that the metric outside the structure tends to FRW all 
the times while the central overdensity region shows a collapse after the initial expansion. 
At the initial conditions, where the density contrast of the overdensity region is still too 
small, we may assume that the metric is almost FRW or LTB; the density contrast and 
the pressure does not play a significant role. The dynamics of Lemaitre universe will give 
us anyhow the expected structure at late times. To choose the initial conditions at the 
time to, we therefore will then use a LTB solution with a negative curvature function. We 
have in fact tried both examples of LTB or FRW initial data and received no significant 
difference between the final Lemaitre solution. 

Now, let us choose the the two initial functions /(r) = f{to,r) and M(r) = M{to,r) in 
the following way to achieve an asymtotically FRW final solution: 

fir) = -^re-^ (18) 

M(r) ^-r^/\l + r^/^). (19) 
a 

Far from the central overdensity region we have 

lirrir^^fir) = 0, (20) 

^3 

limr^^M{r) = — , (21) 
a 

which shows the asymptotically FRW behavior of the initial conditions. The corre- 
sponding LTB solution of Einstein equations now gives us R{r) — Ro{r) at the initial 



time to- Assuming an equation of state is now enough to numerically calculate the 
necessary dynamical functions of the model. Specifically, by looking at R{t, r) and 
p{t, r) we may extract informations of how the central region start collapsing after 
the initial expansion and how a black hole with distinct apparent and dynamical 
horizon develops while the outer region expands as a familiar FRW universe. We 
may also find out the difference to the case of the pressure-less model. It will also 
show us if and how the very weak gravity outside the collapsed structure affects the 
dynamic of the central structure in comparison to the familiar Schwarzschild model. 
The result of the numerical calculation for both equation of states is given in the following. 



A. The density behavior 

The density profile for both equation of states as a function or t and r for the two cases 
of constant w and asymtotically vanishing w is given in Figs.([T]) and (|2]). The development 
of the central black hole and the effect of the pressure on it is evident. Note the more 
significant development of the central pressure in the case of non-vanishing pressure far 
from the center. The over-density region in the collapsing phase is always separated from 
the expanding under-density region through a void not expressible in these figures. We 
will consider the place of the void as the boundary of the structure. This boundary is 
always near to the boundary between the contracting and the expanding region of our 
model. 




FIG. 1: Density evolution of the cosmological black hole for the perfect fluid with the equation 
of state p = wp . 



FIG. 2: Density evolution of the cosmological black hole for the perfect fluid with the equation 
of state p = wps{r). Note the less significant central density and the more flat density profile 
for constant radius near by the center of the structure 

B. The Physical radius 

Figs.(|3]) and (|4]) show the behavior of the velocity R in the more interesting collapsing 
region for the two models. It is evident how the pressure slow down the collapse in the 
central region; With the increase of the pressure the collapse velocity R decreases at the 
center of the black hole. The influence of the pressure on the collapsing velocity near by 
the void is more evident from the behavior of the mass flow described below. 
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FIG. 3: The behavior of the velocity R in the case of p = wp. Evidently the pressure slows 
down the collapse velocity near by the center of the structure. The less the pressure the smaller 
is the radius of the collapsing region 
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FIG. 4: The behavior of the velocity R in the case oi p = wpsir). The features are quahtatively 
as in Fig ([3]). 

C. The Apparent and Event Horizon 

The boundary of a dynamical black hole, where the area law and the black hole tem- 
perature are defined, is a non-trivial concept (see for example jsl and 13|). Our model 
is again a good example to see the behavior of both apparent and event horizon of a 
dynamical structure within an expanding universe. It is easily seen that the apparent 
horizon for our cosmological black hole is located at i? = 2M 10 |. 

This apparent horizon is calculated in t, r coordinates numerically. It is always space- 
like tending to be light-like at late times. This can best be seen by comparing the slope 
of the apparent horizon relative to the light cone at every coordinate point of it. This is 
in contrast to the Schwarzschild black hole horizon where it is always light-like. At the 
late times, however, we expect the apparent horizon to be approximately light-like and 
approaching the event horizon. This is reflected in the Figs.dS]) and ([6]). It is evident 



that 



dr 



AH < 



dt_ 

dr 



\null 



at all times on the apparent horizon, the difference tending to zero 



at late times. Therefore, the apparent horizon is always a space-like dynamical horizon 
leading to a slowly varying horizon at late times . Note that the qualitative result is 
independent of the equation of state. 

We now show how the dynamical horizon of our cosmological black hole becomes a 
slowly evolving horizon at late times. Let's first define the evolution parameter c such 
that the tangent vector to the dynamical horizon, V , is given by 



— cn' 



(22) 



where the two vectors ^°' and are normal null vectors on a space-like two surface S in 
(t, r) plane (see ItI]). We expect c to go to zero at late times in order for our dynamical 



dt 
dr 



null geodesic 




apparent horizon 



FIG. 5: The p = wp case: ^\ah < ^\nuii on the apparent horizon. Therefore, the apparent 
horizon is always a space-like dynamical horizon leading to a slowly varying horizon at late 
times. 
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FIG. 6: The p = ws{r)p case: ^\ah < ^\nuii on the apparent horizon. Qualitatively, there is 
no difference to the other case of the equation of state. 



horizon to become a slowly evolving horizon. In the case of our Lemaitre model c is 
calculated to be 



M' + wM' 



\AH- 



(23) 



M' - wM' - R' ' 

The result of the numerical calculation for different equation of states and different state 
functions is given in Figs. ([7]) and (IH]). We see the decreasing behavior of the function c in 
the course of time independent of the equation of state. We may then conclude that the 
dynamical horizon of the cosmological black hole tends to a slowly evolving horizon. 




FIG. 7: The p = wp case: the more pressure the sooner the dynamical horizon becomes a slowly 
evolving horizon. 




FIG. 8: The p = wps{r) case: qualitatively, the same behavior as the other case. 



D. Mass and matter flux 



Due to the expanding background we expect the matter flux into the djTiamical black 
hole to be decreasing and the dynamical horizon to become a slowly evolving horizon lol |. 
We know already that there is no unique concept of mass in general relativity correspond- 
ing to the Newtonian concept. 

The question of what does general relativity tell us about the mass of a cosmological 
structure in a dynamical setting was discussed recently jl^. It was shown H] that The 
Misner-Sharp quasi-local mass, M, is very close to the Newtonian mass. 

Let us then take the Misner-Sharp mass for this black hole and calculate the corre- 
sponding matter flux into the black hole. In the case of Lemaitre model, the matter flux 



according to the comoving time is given by 
dM{r, t) 



dt '^^ dt ' dr dt 

The result of the numerical calculation is depicted in Figs. ([9]) and (fTOj) . Note how the 
pressure decreases the rate of matter flux into the black hole. 



9M(r, t) I (9M(r, t) dr 

AH H 7{ 7rr\AH 



dv 

M\ah + M'—\ah (24) 
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FIG. 9: The p = wp case: the rate of matter flux into the black hole decreases with the pressure. 
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FIG. 10: The p = ws{r)p case: qualitatively the same behavior as the other case. 



IV. DISCUSSION 

We have studied the evolution of a structure made of perfect fluid as an exact solution of 
Einstein equations within an otherwise expanding FRW universe. The structure boundary 
is marked through a void from the expanding part of the model which is very much close 
to a FRW universe already near by the void. The collapsed region develops to a dynamical 
black hole with an apparent horizon tending to a slowly evolving horizon at late times, 



where the mater flux into the black hole also decreases preventing more matter falling 
inside the black hole. The space-like apparent horizon tends approximately to a light-like 
event horizon at late times. We have, therefore, to conclude that the mere existence of 
a cosmological matter, even dust, may have significant effect on the central black hole 
differentiating it from a Schwarzschild one irrespective of how small the density outside the 
structure is. Hence we may not be allowed to speak about the Newtonian approximation 
because of the very weak gravity in cases of non-local quantities such as the horizon and 
the mass. 
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